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Abstract. TeV 7-ray emission has been recently discovered by Cherenkov telescopes from two microquasars, 
LS 5039 and LS I +61°303. This emission is likely to be produced inside the binary system since in both cases 
variability with the orbital period of the binary has been discovered. In fact, such emission features have been 
recently predicted by the inverse Compton (IC) pair cascade model. In this model, electrons accelerated in 
the jet develop a cascade in the anisotropic radiation of the massive star. The 7-ray spectra emerging from the 
cascade strongly depend on the location of the observer with respect to the orbital plane of the binary. Here 
we apply this model in order to investigate the possible 7-ray emission features from another compact massive 
binary of the microquasar type, Cyg X-l. We conclude that the observational constraints at lower energies (from 
MeV to GeV) suggest that the spectrum of electrons injected in the jet is likely steeper than in other TeV 7-ray 
microquasars. The cascade 7-ray spectrum produced by electrons with such spectrum in Cyg X-l should be below 
the sensitivities of the MAGIC and VERITAS class Cherenkov telescopes. However, if the electron spectrum is 
flatter, then the highest TeV 7-ray fluxes are predicted at ~7 hrs before and after the phase when the compact 
object is in front of the massive star. We suggest that Cherenkov telescopes should concentrate on these range of 
phases since the TeV flux can vary by a factor of ~20 with the period of the Cyg X-l binary system. Moreover, 
the model predicts clear anticorrelation of the GeV and TeV 7-ray emission. This feature can be tested by the 
future multiwavelength observations with the AGILE and GLAST telescopes in the GeV energy range and the 
MAGIC and VERITAS telescopes in the TeV energy range. 
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1. Introduction 

The importance of high energy processes in X-ray binary 
systems has been hypothesized for about two decades and 
in fact some evidence (although not statistically signifi- 
cant) has been reported in the literature (for review see 
e.g. Weekes 1988, 1992). However, only one of these early 
Cherenkov observations of Cen X-3 has been later con- 
firmed, and remains unquestioned up to the present time 
(Chadwick et al. 2000, Atoyan et al. 2002). Some of these 
well known binary systems have been also related to the 
EGRET sources above 100 MeV due to directional coin- 
cidence since they are close or inside of relatively large 
error boxes of some EGRET sources: LS I +61 303 = 
2CG135+01 = 2EG J0241+6119 (Thompson et al., 1995, 
Kniffen et al. 1997), Cyg X-3 = 2EG J2033-4112 (Mori et 
al. 1997), Cen X-3 (Vestrand, Sreekumar & Mori, 1997), 
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and LS 5039 = 3EG J1824-1514 (Paredes et al. 2000). 
However, only very recently some binary systems belong- 
ing to the class of high mass binaries, have been explic- 
itly detected at high energy 7-rays by Cherenkov tele- 
scopes. The HESS collaboration reported the detection 
of TeV 7-ray emission from the Be type binary system 
PSR B1259-63, close to the periastron passage (Aharonian 
et al. 2005). Also TeV 7-rays have been detected from 
LS 5039 (HESS collaboration, Aharonian et al. 2005, 
Aharonian et al. 2006) and LS I +61 303 (MAGIC col- 
laboration, Albert et al. 2006). These two objects belong 
to the class of the so-called microquasars. These objects 
show evidences of anisotropic outflows usually interpreted 
as jet structures similar to the rclativistic jets observed in 
active galactic nuclei. Such jets in microqusars might be 
created in the inner parts of the accretion disks around 
compact objects such as a neutron or a solar mass black 
hole that is accreting matter from the companion star. 
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In fact, microquasar features, very similar to those ob- 
served in LS 5039 and LS I +61 303, are also observed in 
the well-known object, Cyg X-l. It has been detected by 
the COMPTEL instrument up to ~ 10 MeV with a spec- 
trum very close to a power law with differential spectral in- 
dex 2.6 (McConncll et al. 2000, 2002). However, the source 
has not been observed at higher energies by the EGRET 
instrument (sec the upper limit in McConncll et al. 2000). 
The MeV 7-ray emission is usually interpreted in terms of 
two component thermal and non-thermal comptonization 
models from the accretion disk (e.g. Gicrlihski et al. 1999), 
or as a result of comptonization of the stellar photons by 
relativistic electrons inside the jet (Bednarek et al. 1990, 
Georganopoulos et al. 2002, Romero et al. 2002). The gen- 
eral picture of the production of 7-rays is that they arise 
as the result of comptonization of stellar photons by elec- 
trons in the jet. This model has been recently applied to 
the microquasar LS 5039 (Dormer & Bottcher 2006 and 
Bosch- Ramon, Romero & Parcdes 2006). However, these 
models do not take into account Inverse Compton (IC) e 
pair cascade processes. Cascade processes become impor- 
tant for these binary systems since the optical depth for 
electrons (which are isotropic in the jet) depends largely 
on their temporary directions (for detailed calculations 
of the optical depths, see Bednarek 1997, 2000 and more 
recently e.g. Sierpowska & Bednarek 2005, Bottcher & 
Dermer 2005, Dubus 2006, Bednarek 2006a). Due to this 
effect, electrons lose energy mainly in directions for which 
the optical depths are the largest, but not in the direc- 
tion toward the observer (see the previously cited pa- 
pers). The calculations of the 7-ray spectra, which are 
produced in the anisotropic IC pair cascade toward 
the observer, and their light curves as a function of the 
phase of the binary system, have been recently calcu- 
lated by (Bednarek 2006a). These calculations are based 
on the previous work done for other compact binary sys- 
tems such as Cyg X-3 (Bednarek 1997), and Cen X-3 
(Bednarek 2000, Bednarek & Giovannelli 1999). Recently, 
the cascade model for the high energy processes inside 
compact massive binaries has been successfully applied to 
the binary system LS I +61 303, predicting the phases at 
which the maximum TeV 7-ray emission can be expected 
(Bednarek 2006a, b). Such a prediction was confirmed by 
MAGIC observations (see Albert et al. 2006). As an alter- 
native scenario, hadronic production of 7-rays and neu- 
trinos from these binaries has been also considered (e.g. 
Romero, Christiansen & Orellana 2005, and Aharonian et 
al. 2006). 

In view of the recent discoveries of the TeV 7-ray emis- 
sion from microquasars LS 5039 and LS I +61 303, we 
consider the IC pair cascade model for possible high 
energy processes in another binary system, Cyg X-l/HDE 
226868. 

2. The binary system Cyg X-l/HDE 226868 

Cygnus X-l has been the most extensively studied stel- 
lar mass black hole candidate since its discovery as one 



Table 1. Parameters of the binary system and the jet. 



Parameter 


Value 


R+: stellar radius 


20 R Q 


T, : surface temperature 


3 x 10 4 K 


M^. stellar mass 


40 ± 10 M Q 


M x : companion mass 


16 + 3 Mq 


a: orbital radius 


2.157?, 


e: eccentricity 





i: inclination 


30° 


rj n : disk inner radius, radius of the jet 


10 ' cm 


a: jet opening angle 


0.1 rad 


jZmin: inner boundary of the active jet 


0.1R-, 


z max : outer boundary of the active jet 


WR* 


r\: magnetization parameter of the jet 


1.0-0.01 


Bd'- magnetic field at the base of the jet (for r\ = 1) 


10 s G 


£: acceleration efficiency of electrons 


0.5-0.01 



of the most intense X-ray sources in the sky (Bowyer et 
al. 1965). In the X-rays the source shows two emission 
states (see, e.g. Pricdhorsky, Terrell & Holt 1983, Ling et 
al. 1983, Liang & Nolan 1983), spending most of the time 
(~90%) in the so-called hard state (low level of soft X-ray 
emission) and the remaining of the time in the soft state. 
Cyg X-l was also identified as a radio source by Braes 
& Miley (1971). High resolution radio observations show 
a steady jet from this source in the hard X-ray spectral 
state (Stirling et al. 2001) and also a radio and optical 
lobe oriented along the jet direction (Gallo et al. 2005). 
The jet appears when the X-ray spectrum starts softening 
and its intrinsic velocity is estimated on > 0.3c (Fender et 
al. 2006). 

Hjellming & Wade (1971) identified Cyg X-l with the 
optical star HDE 226868, which is an 09.7 lab type mas- 
sive star inside a binary system with a period of 5.6 days 
(Webster & Murdin 1972, Bolton 1972). The main param- 
eters of this binary system have been recently confirmed 
by Ziolkowski (2005): the mass of the optical companion 
M* = 40 ± 10 M Q ; its radius R* = 22 ± 3 R Q ; an effective 
temperature of T, = 30, 000 ± 2, 000 K for a distance to 
the system D = 1.8-2.35 Kpc; and a mass of the compact 
object M x = 16+3 M Q which is in an almost circular orbit 
e = 0; with an inclination angle of the binary i = 33° ±5°. 
In this paper we adopt the values for the binary system 
shown in Table 1. 

3. A microquasar in the compact massive binary 

We consider a generally accepted scenario for micro- 
quasars, in which electrons are accelerated continuously 
along the jet. The jet propagates from the compact ob- 
ject along an axis perpendicular to the plane of the bi- 
nary system. The details of such a general scenario have 
already been described in Bednarek (2006b). Here we re- 
peat the basic points of the scenario for the sake of com- 
pleteness. The shock acceleration mechanism is usually 
accepted as being responsible for the acceleration of the 
electrons. It is assumed that electrons are injected with 
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the power law spectrum extending up to the maximum 
energies determined by the balance between the accelera- 
tion rate and the most efficient energy loss rate. In the case 
of microquasars, the synchrotron process (determined by 
the magnetic field strength locally in the jet) or the inverse 
Compton process (determined by the radiation field of the 
massive star) usually provide the most efficient cooling for 
the electrons. 

The magnetic field strength inside the jet can be es- 
timated on the basis of the cquipartition arguments con- 
cerning the relation between the accretion disk luminos- 
ity and the kinetic power of the jet, i.e. the so called 
jet-accretion-disk symbiosis model developed by Falcke et 
al. (1995), Falcke et al. (2004). In fact, in the case of Cyg 
X-l, the X-ray luminosity (probably produced in the ac- 
cretion disk) and the jet kinetic power are estimated to be 
of the same order (Gallo et al. 2005). Note, however, that 
this is not always the case since some other binaries are 
clearly undcrluminous, e.g. LS 5039 and LS I +61 303. 
Moreover, it is commonly expected, on the basis of the 
cquipartition arguments, that the energy density of the 
magnetic field inside the jet should be related to the ki- 
netic power of the jet (see e.g. Bosch Ramon et al. 2005). 
Therefore, we accept the hypothesis that the magnetic 
field energy density in the inner part of the accretion disk 
(at the base of the jet) should be related to the energy 
density of the disk radiation (see also e.g. Bednarek 2005). 
On the basis of these arguments we arrive to the conclu- 
sion that the magnetic field in the jet can be as large as 
£?d ~ 10 5 G (e.g. Bednarek 2006b). Assuming a simple 
conical structure for the jet (fixed opening angle a), the 
magnetic field in the jet should drop with the distance 
from its base according to 

B(x) = B d T)/(l + az) rj B d ?]/(az) for z > 1/a, (1) 

where a is assumed of the order of 0.1 rad, rj is the ratio 
of the energy density of the magnetic field to the energy 
density of the disk radiation at the base of the jet (mag- 
netization parameter of the jet). This parameter describes 
how far the magnetic field is from the cquipartition with 
respect to the kinetic energy density of the jet plasma, z is 
the distance along the jet from its base expressed in units 
of the disk inner radius n n , where r; n is of the order of 10 7 
cm for the compact object with the mass of the order of 
10 M Q . 

Having defined the magnetic field strength at different 
parts of the jet, we can estimate the acceleration rate of 
the electrons (with energy E and the Lorentz factors 7), 

^acc(7) = &E/r L rj 10 13 6B eVs-\ (2) 

where £ is the acceleration efficiency, = E/eB is their 
Larmor radius, B is the magnetic field at the accelera- 
tion site (in Gauss), e is the electron charge, and c is the 
velocity of light. In Eq. 1, a fundamental role is played 
by the acceleration efficiency of particles £, but this is in 
fact unknown. There are some observational suggestions 
that £ should be not far from unity in the case of rel- 
ativistic shocks, e.g. presence of relativistic leptons with 



energies ~ 10 15 eV inside the Crab Nebula. On the other 
hand, the theory of shock acceleration processes estimates 
£ as - 0.1(w/c) 2 , i.e. > 0.01 for v > 0.3c (Malkov & 
Drury 2001). Therefore, in this paper we consider the val- 
ues of the acceleration efficiency of electrons in the jet 
in the range £ = 0.5 — 0.01. The basic parameters which 
describe the jet are shown in Table 1. 

As we noted above, during the acceleration process, 
electrons lose energy mainly by the synchrotron and the 
inverse Compton (IC) processes. We are interested in sit- 
uations in which electrons can be accelerated to TeV en- 
ergies. This is not easy to obtain relatively close to the 
base of the jet where the magnetic field is large and the 
synchrotron losses very strong. Therefore, since we are in- 
terested in efficient injection of high energy IC 7-rays by 
electrons in the jet, it makes sense to consider distances 
from the accretion disk for which the IC energy losses on 
the massive star radiation are at least comparable to syn- 
chrotron losses. In this cascade model we will concentrate 
on the parameter range adequate for the production of 
TeV 7-rays. 

Let us estimate the maximum energies of electrons 
in the shock acceleration process allowed by the energy 
losses. In these estimations we concentrate on the parts 
of the jet which are above z m [ n = O.IRq = 1.4 x 10 11 
cm, since very close to the base of the jet the radia- 
tion field from the accretion disk (and disk corona) may 
dominate over the massive star radiation field (total lu- 
minosity of the disk is estimated in several times 10 37 
erg s _1 , see e.g. Poutancn, Krolik & Ryde 1997). At 
z m i n , the disk radiation field is diluted by the factor of 
~ (/3z min /r in ) 2 rj 10~ 6 - 10~ 4 , where r in ~ 10 7 cm is 
the inner radius of the disk around the 10M Q black hole, 
and the parameter (3 = 7"i n /r cor ~ 0.01 — 0.1 defines 
the extend of the disk corona assumed in the range of 
r cor rj 10 - 100rj n . Note that the massive star luminos- 
ity in Cyg X-l is clearly larger than the disk luminosity. 
Therefore, disk radiation can become important only in a 
relatively small part of the inner jet. 

The energy loss rate by the synchrotron process is 
given by, 

A y „(7) = ^vtcUbJ 2 « 2.7 x HT 3 B 2 7 2 eVs" 1 , (3) 

where Ub = B 2 /8ir is the energy density of the magnetic 
field, 7 is the Lorentz factor of the electrons, and ctt is 
the Thomson cross section. Electrons lose energy on the 
IC process in the Thomson (T) regime, if their Lorentz 
factors are, 

7 « 7x/kN - mc 2 /3k B T «2x 10 5 /r 4 , (4) 

(where T = 10 4 T4 K is the surface temperature of the 
massive star which defines the black body spectrum, and 
/cb is the Boltzman constant), and in the Klein-Nishina 
(KN) regime, for the Lorentz factors 7 ^> 7t/kn- The 
energy loss rate in the Thomson regime is given by 

Ac(7) = ^^Tci7 r ad7 2 «3.8T 4 4 7 2 A 2 eVs" 1 , (5) 
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log z log z log z 

Fig. 1. The maximum Lorentz factors of electrons as a function of the distance, z from the base of the jet (measured 
in units of the stellar radius R±) for different acceleration efficiencies: £ = 0.5 (from the upper, thickest curve), 0.1, 
0.03, and 0.01 (to the bottom, finest curve) and two ratios of the magnetic field to disk radiation energy densities at 
the base of the jet equal to ?/ = 1 (a), 0.1 (b), and 0.01 (c). The Maximum Lorentz factors are obtained from the 
comparison of the acceleration rate with the energy loss rate on synchrotron process (dotted lines), inverse Compton 
process on scattering of stellar radiation from the massive star HDE 226868 in the Thomson regime (full curve) and 
the Klein-Nishina regime (dashed curves). The base of the jet is at the distance of 2.15.R* from the center of the 
massive star. 



where t/ rad = 4.7 x 10 13 T 4 4 /r 2 eV cm" 3 , and r is the dis- 
tance to the center of the massive star in units of its radius 
R*. The energy loss rate in the KN regime depends only 
logarithmically on the Lorentz factor of the electrons. We 
approximate these losses by (see e.g. Blumenthal & Gould 
1970), 

Ac N (7) « Ac(7t/kn) In (4fc B r 7 /mc 2 - 2). (6) 

The maximum energies of the accelerated electrons are de- 
termined by the balance between the acceleration mecha- 
nism and by the most efficient mechanism for the energy 
losses. The comparison of Eqs. [2] and [3] gives the maxi- 
mum allowed Lorentz factors of the electrons due to the 
synchrotron energy losses, 

7sy „ « 6 x lO^/fl) 1 ' 2 . (7) 

The maximum energies allowed by the IC process in the 
Thomson regime are (from Eqs. [2] and [5]), 

7l T c « 1.6 x IQ^iBf^r/Tl (8) 

and in the KN regime arc (from Eqs. [2] and [6]), 

7l K c N a 1.5 x 10 5 (2 + e^ Br2 ^)/T 4 . (9) 

Note that the limit given by Eq. [8] is only valid provided 
that < 7t/kn; i- c when the parameters of the jet and 
the massive star fulfill the following condition, T^/r > 
8(£-B) ' 5 (from Eqs. [Hand [8]). In any other case, the limit 
given by Eq. [5] has to be taken. 

The Larmor radii of electrons accelerated to the maxi- 
mum Lorentz factors allowed by the radiation constraints 
has to be smaller than expected from the constraints im- 
posed by the size of the jet. The Larmor radius of elec- 
trons (defined below Eq. 2) is thus lower than the size 
of the jet, tl < (1 + az)r- m (see above), provided that 



7size < eBdTi/m c . For the considered parameters and the 
applied model for the magnetic field distribution along the 
jet (see Eq. 1), 7 S ; ZC < 6 x 10 8 . Note that this limit is in- 
dependent of the distance from the base of the jet. It does 
not impose additional constraints on the values obtained 
above for distances from the base of the jet which are of 
interest in this model. 

The maximum Lorentz factors to which electrons can 
be accelerated at a specific part of the jet are deter- 
mined by the acceleration efficiency and the energy losses 
through synchrotron and IC processes. These energy loss 
processes depend on the conditions in the jet (i.e., mag- 
netic field strength and acceleration efficiency) , on the pa- 
rameters of the massive star (its radius and surface tem- 
perature), and on the distance of the acceleration region 
from the massive star. In Fig. [TJ we plot the maximum 
Lorentz factors of electrons as a function of the distance 
from the base of the jet (given by either Eq. [71 or Eq. 
or Eq. [9]) , for selected values of the acceleration efficiency, 
£, and the ratio, rj. The maximum energies of the acceler- 
ated electrons increase with the distance from the base of 
the jet as oc z 1 / 2 , when the acceleration process is satu- 
rated by synchrotron energy losses (see also Bosch Ramon 
et al. 2006). However, in the regions of the jet where the 
magnetic field is relatively weak (farther from the base of 
the jet but still close to the massive star), the accelera- 
tion process is saturated by IC energy losses. Then, the 
maximum Lorentz factors of the electrons are significantly 
lower (note the broad parabolic dips in Fig. 1, especially 
evident during the periastron passage of the compact ob- 
ject). 

In a previous work (Bednarek 2006b), we have shown 
that electrons accelerated in the jet have to cool locally, i.e. 
close the the acceleration site, due to the very large energy 
losses, provided that they are injected within ~ 100i?* 
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Fig. 2. (a) The optical depths for 7-rays (as a function of their energy) on pair production in collisions with stellar 
photons calculated for 7-rays propagating from a specific distance from the massive star up to the infinity. 7-rays are 
injected at three different distances from the massive star HDE 226868, z = 2R+ (a), z = 5R+ (b), and z = 20i?* (c), 
and for selected angles measured from the direction defined by the centers of the stars, a = 0° (i.e. toward the center 
of the massive star the optical depth is calculated only to the surface of the star, thick full curve), 30° (thick dashed), 
60° (thick dotted), a = 90° (thick dot-dashed), 120° (thin full), 150° (thin dashed), and 180° (thin dotted). 



from the base of the jet. At larger distances accelerated 
electrons are advected outside the binary system (along 
the jet) and can produce high energy radiation through 
other processes. However in such a case, we do not expect 
any clear variability for the 7-ray emission with the period 
of the binary system. Therefore, we do not consider that 
possibility in this paper. 

Following our previous work, we also assume that 
electrons are injected with the power law spectrum de- 
scribed by the single spectral index independent on the 
distance from the base of the jet z (typical values from 
the range 2-2.6 arc considered). However, the injection 
rate of electrons, i.e. the power transfcrcd from the mag- 
netized plasma to rclativistic electrons L c , can differ, de- 
pending on a range of possible processes occurring in the 
jet. Two different likely scenarios are considered: 

1. constant injection of electrons along the jet from its 
base at z mm : L c = const., 

2. injection rate drops along the jet according to L c (z) oc 
z-\ 

3.1. Propagation of 7- rays 

Since the Cyg X-l binary system is very compact (with 
a separation of the stars of only a = 2.15i?*), 7-rays pro- 
duced in the part of the jet which is closer than a few tens 
of R + from the massive star can be efficiently absorbed 
by collisions with the thermal radiation from the massive 
star. Following our previous works (Bednarek 1997, 2000, 
Sicrpowska & Bednarek 2005), we calculate the optical 
depths for 7-rays as a function of the energy in the most 
general case, i.e. for an arbitrary distance from the mas- 
sive star and direction of their propagation (for an angle 
measured from the direction defined by the centers of the 
stars). The results are shown in Fig. [2] for an injection 
distance of 7-rays from the center of the massive star at 



2i?* (a), 5i?+ (b), and 20-R* (c). It is clear that the opti- 
cal depths are significantly above unity at specific direc- 
tions even if the injection distance from the massive star 
is quite large. Since the cross sections for the IC process 
in the KN regime and for the 7+7 collisions (7+7 — * e ± ) 
are comparable, it is expected that relativistic electrons 
injected isotropically inside the jet lose energy mainly in 
directions where the optical depths are the largest. They 
inject primary 7-rays mainly in directions of the largest 
optical depths, i.e. toward the source of soft radiation. 
These primary 7-rays can develop IC e pair cascades 
in the anisotropic radiation field of the massive star. A 
part of these cascade 7-rays escape toward the observer 
located at specific directions with respect to the binary 
system. We follow such a complicated cascade process ap- 
plying the Monte Carlo cascade code described in detail in 
our previous works (e.g. Bednarek 2000, 2006a). The code 
includes also the effects connected with synchrotron cool- 
ing of the electrons injected into the magnetic field of the 
jet. Therefore, we are also able to calculate the simulta- 
neous synchrotron photon spectra produced by electrons 
accelerated in the jet. 7-ray photons escaping from the bi- 
nary system are sorted in specific ranges of solid angles 
on the hemisphere and grouped at specific ranges of en- 
ergies. Therefore, as a final product, we obtain the 7-ray 
spectra produced in such anisotropic cascades at specific 
directions on the hemisphere. 

We note that in parts of the jet relatively close to the 
accretion disk (this part of the jet is not considered in the 
present paper), 7-rays can also be absorbed in the disk 
radiation field (see, e.g., the calculations of the optical 
depths by e.g. Herterich 1974, Protheroe & Stanev 1987, 
Carraminana 1992, Bednarek 1993). The disk radiation 
can be also upscattered to the 7-ray energies by electrons 
provided that they are accelerated very close to the base 
of the jet as considered by e.g. Aharonian et al. (1985) or 
Romero et al. (2002). 
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Fig. 3. Cascade 7-ray spectra produced by electrons injected at different distances from the base of the jet z = O.lr* 
(full curve), li?* (dashed), and lOr* (dotted) for two locations of the compact object on its orbit around the massive 
star, in front of the massive star (upper figures) and behind the massive star (bottom figures). The 7-ray spectra 
expected from the jet moving inside the hemisphere containing the observer, who is located at the inclination angle 
of i = 30°, are shown on the left figures, from the counter-jet (middle figures), and the sum of the jet + counter-jet 
spectra (right figures). The acceleration of electrons in the jet is defined by the parameters £ = 0.1 and 77 = 0.1 and 
the electrons are injected with a power which drops along the jet according to L c (z) oc z~ 2 . 



4. Cascade gamma-ray spectra 

We investigate the features of the 7-ray spectra produced 
in such cascade process as a function of parameters de- 
scribing the acceleration of electrons in the jet. Since the 
Doppler factor of the jet is unknown in the case of Cyg 
X-l (the lower limit on the jet velocity is ~ 0.3c - Fender 
et al. (2006) - and the inclination angle of the binary sys- 
tem is i ~ 30°), two cases should be considered. When 
the Doppler factor is low, as e.g. in the case of LS I +61 
303 (Bednarek 2006b), then we have to consider the 7- 
ray spectra produced not only by electrons in the jet (di- 
rected toward the observer), but also the 7-rays produced 
by electrons in the counter-jet (directed in the opposite 
direction). However, for the relativistic jets with Doppler 
factor D 3> 1, it is enough to consider only 7-ray emis- 
sion from the observer-directed jet since the contribution 
from the counter-jet is suppressed due to relativistic ef- 
fects. Below we study the basic features of the cascade 
7-ray spectra, i.e. their dependence on the distance from 
the base of the jet, on specific parameters defining the ac- 
celeration of electrons, and their dependence on the phase 
of the binary system and the spectra at specific phases. 
In the next section we discuss the GeV and TeV 7-ray 
spectra which might be expected from the massive binary 
system Cyg X-l. 



4.1. Dependence on the distance from the base of the 
jet 

Let us at first calculate the 7-ray spectra escaping from 
the binary system toward an observer located at the in- 
clination angle i = 30°, These 7-rays are produced by 
electrons at different parts of the jet, and their produc- 
tion is parameterized by the distance from the base of the 
jet. The parameters defining the acceleration of electrons 
have been chosen from the center of the range of values 
specified above (i.e. £ = 0.1 and 77 = 0.1 - see Fig. [T]). 
For easier interpretation of the results, we apply the dif- 
ferential power law spectrum of the injected electrons of 
the type oc E~ 2 (equal power per decade) up to the max- 
imum energies defined by £ and 77. In Fig. [31 we show 
the cascade-produced 7-ray spectra (multiplied by the en- 
ergy squared) that are produced by electrons in the jet, 
the counter-jet, and the sum of both of them, i.e. jet + 
counter-jet spectra. The upper and lower panels show the 
spectra which are produced when the compact object is 
in front of or behind the massive star, respectively. The 
shape of the 7-ray spectra clearly depend on the injection 
distance of electrons. 

Electrons which are injected closer to the base of the 
jet are also closer to the surface of the massive star and 
therefore are immersed in a stronger soft radiation field. 
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These cascade 7-ray spectra show a low level of TeV emis- 
sion, due to the absorption of 7-rays, and a relatively high 
level of GeV emission (produced due to the efficient cas- 
cading). Significant fluxes of TeV 7-rays can only be pro- 
duced in the upper part of the considered jet. The cut-offs 
in the 7-ray spectra are determined by the maximum en- 
ergies of electrons accelerated in the jet. These energies 
increase with the distance from the base of the jet. The 
7-ray spectra produced toward the observer in the cascade 
process initiated by electrons in the counter-jet typically 
show higher GeV but lower TeV emission. This is due 
to the fact that these 7-rays arrive to the observer after 
suffering from a severe cascade process (the angle taken 
from the direction between the jet-counter-jet axis and 
the direction toward the observer is much larger). In fact, 
electrons in the counter-jet barely contribute to the 7-ray 
spectrum above ~ 1 TeV. Therefore, the sum of the jet + 
counter-jet spectra (right figures) is similar to the spectra 
from the jet with some increase in the GeV energy range 
due to the counter-jet contribution. The counter-jet only 
has the effect of changing the normalization of the total 
escaping spectra in the GeV energy range. 

4.2. Gamma-ray spectra integrated along the jet 

In Fig. IU we investigate the dependence of the cascade 
7-ray spectra escaping toward the observer as a function 
of the parameters defining the acceleration model (see the 
profiles of the maximum electron energies in Fig. [1} . The 
7-ray spectra produced in terms of this model are inte- 
grated over the part of the jet extending between (0.1- 
10)R^, since only from this limited part of the jet can the 
cascade emission dominate (in the radiation of the mas- 
sive star) and the 7-ray emission be correlated with the 
orbital period of the binary system (as recently discovered 
in the case of LS I +61 303, Albert et al. 2006). 

At distances closer than 0.1R*, the energy losses via 
synchrotron process and scattering of the disk radiation 
might dominate if all electrons can be accelerated in that 
region to above ~1 TeV. At distances greater than 10R*, 
electrons accelerated in the jet are not able to cool locally 
in the jet. They are advected at large distances along the 
jet due to the motion of the plasma inside the jet, where 
they cool via other processes. Note, also, that scenarios in 
which the possible production of 7-rays occurs very close 
or far away from the base of the jet do not provide an 
explanation of the correlation of the 7-ray emission with 
the period of the binary system. 

From a comparison of calculations performed for differ- 
ent parameters, we conclude that strong TeV 7-ray emis- 
sion appears for larger values of the acceleration efficiency 
of electrons in the jet (full and dotted curves in Fig. [4|). 
As expected, the 7-ray spectra generally cut off at lower 
energies for the case of less magnetized jets and lower ac- 
celeration parameters. This is due to the saturation of 
electron acceleration at lower energies by IC energy losses 
in these cases (see Fig. [1]). The 7-ray spectra produced 
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Fig. 4. The cascade 7-ray spectra (multiplied by energy 
square) produced by electrons in the jet (left figures) and 
in the jet + counter-jet (right figures), for the situations 
when the compact object is behind (upper figures) and in 
front of the massive star (bottom figures). Specific spectra 
have been obtained for: £ = 0.1 and 77 = 0.1 (full curves), 
£ = 0.01 and n = 0.1 (dashed), £ = 0.1 and n = 0.01 
(dot-dashed), £ = 0.3, and r\ = 0.1 (dotted). They are 
integrated over the part of the jet extending between 
z=0.1-10R* measured from the base of the jet. Electrons 
are injected with the power law spectrum, oc E~ 2 . The 
power in electron spectrum drops along the jet according 
to LJz) oc z~ 2 . 



in the cascade initiated by electrons, when the jet is be- 
hind the massive star, show a sharp cut off at ~ 100 GeV 
due to the efficient cascading for all the considered cases. 
However, when the compact object is in front of the mas- 
sive star, the TeV 7-ray emission can reach the observer 
with a relatively flat spectrum (spectral index ~ 2.0-2.5 
above a few 100 GeV) and the flux at a level of ~0.1-0.01 
of the GeV emission. We conclude that the 7-ray emission 
produced by electrons in the jet in terms of the considered 
model has very interesting features which should allow us 
to explicitly distinguish between different models for the 
7-ray production in massive binaries, e.g. scattering of the 
disk radiation by electrons in the jet (as considered in the 
case of active galaxies by, e.g. Dernier & Schlickciscr 1993) 
or production of 7-rays far away from the base of the jet 
(e.g. Atoyan & Aharonian 1999). 

A3. Gamma-ray light curve and phase dependent 
spectra 

Let us now discuss the dependence of 7-ray emission as a 
function of the phase of the binary system, where the zero 
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Fig. 5. Gamma-ray light curves in the energy range 1-10 GeV (dashed curves) and >200 GeV (full curves) produced 
by electrons accelerated in the jets of Cyg X-l for the model defined by the parameters £ = 0.1 and rj = 0.1, the power 
law spectrum of electrons, oc E~ 2 , and the power transferred to electrons from the acceleration mechanism which 
drops along the jet according to L c (z) oc z~ 2 . 7-rays produced by electrons in the jet (thin curves) and in the jet + 
counter-jet (thick curves). The middle and right panels show the cascade IC 7-ray spectra escaping to the observer 
from the binary system for three phases of the compact object (measured from the location of the compact object 
behind the massive star): 0. (dashed curve), 0.45 (dotted), and 0.5 (full). 



phase is taken when the location of the compact object is 
behind the massive star. Since the orbit of the compact 
object is almost circular, the 7-ray light curves should be 
symmetrical. We calculate the number of 7-ray photons 
that escape from the binary system toward the observer 
located at the inclination angle i = 30°, in two energy 
ranges: the so-called GeV energy range (between 1-10 
GeV available for the past EGRET telescope and the fu- 
ture AGILE and GLAST telescopes), and the TeV energy 
range (>200 GeV) available for the present and future 
Cherenkov telescopes on the northern hemisphere such as 
the MAGIC and VERITAS telescopes. As an example, we 
concentrate on the model defined by £ = 0.1 and 77 = 0.1 
(the parameters from the center of considered range) , and 
the electron spectrum described by oc E~ 2 , and its power 
dropping along the jet according to L e (z) oc z~ 2 . The 
7-ray light curves (number of photons arriving to the ob- 
server in the GeV and TeV energy ranges) are shown in 
Fig. [5] for the case of 7-ray emission only from the jet 
(thin curves) and the sum of the emission from the jet 
+ counter-jet (thick curves). The TeV 7-ray emission can 
change in this case by more than a factor of 20, showing 
the minimum when the compact object is behind the mas- 
sive star and a sharp maximum for the phases at ^0.05 
before and after the situation when the compact object is 
in front of the massive star. 

Exactly in front of the massive star, the TeV 7-ray flux 
becomes a factor of ~2 lower. Due to the lowest optical 
depths for 7-rays and leptons production of 7-rays is less 
favourite in such a geometry than for surrounding phases. 
The contribution of the counter-jet to the maximum in the 
TeV 7-ray light curve is at the level of ~ 30% of the jet 
emission (compare the thin and thick full curves). On the 
other hand, the GeV 7-ray light curve behaves differently. 
It has the maximum when the compact object is behind 
the massive star and minimum exactly when the compact 
object is in front of the massive star. The maximum is due 



to the fact that the 7-ray spectrum is formed in the very 
efficient cascade process when the compact object and the 
observer are on the opposite sites of the massive star. The 
contributions of the jet and the counter-jet to the GeV 7- 
ray light curve are comparable. It is remarkable that there 
is a very interesting anticorrelation between the 7-ray light 
curves in the GeV and TeV energy ranges. This is the 
unique feature which can serve as an important diagnostic 
of this model for the high energy processes occurring in 
the 7-ray emitting massive binary systems. 

Fig. [3] also shows the 7-ray spectra emerging toward 
the observer for some characteristic phases of the binary 
system for the case of the acceleration of electrons in the 
jet and the sum for the jet + counter-jet emission. Note 
the dependence between the spectral shape and the flux 
of the TeV 7-ray emission. It is predicted that the lower 
level of 7-ray flux should be connected with a flattening of 
the spectrum above ~200 GeV and vice versa. The contri- 
bution of 7-ray emission from the counter-jet has an effect 
of increasing the GeV 7-ray flux with no effect on the TeV 
flux. However, due to the higher GeV emission, the rela- 
tive power in the GeV and TeV energy ranges changes. 
This has consequences for predicting the TeV flux from 
the normalization to the GeV (EGRET) observations. 

5. The multifrequency spectrum of Cyg X-l 

The 7-ray spectrum of Cyg X-l is in fact very poorly con- 
strained. The spectrum in the soft X-ray state shows the 
power law tail with the differential spectral index -2.6 ex- 
tending up to ~10 MeV, without evidence of a cut off 
(McConnell et al. 2000, 2002). However, the jet appears 
only in the low-hard state, so this emission can be only 
considered as an upper limit on the MeV 7-ray flux re- 
lated to the jet. The upper limit at higher energies (> 100 
MeV) has been derived (see McConnell et al. 2000), based 
on the 3rd EGRET catalog (Hartman et al. 1999). 
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Fig. 6. The multiwavelength spectrum of Cyg X-l (differential spectrum multiplied by energy squared) is compared 
with the calculated cascade 7-ray spectra for three acceleration models of electrons in the jet: (1) differential spectral 
index of electrons oc E~ a and a = 2, and the power transfered to electrons changes with distance along the jet 
according to L c (z) oc z~ 2 (model A); (2) a = 2.6 and L e (z) oc z~ 2 (model B), (3) a = 2.6 and L c (z) = const (model 
C). The cut-off in the electron spectrum is determined by £ = 0.1, J] = 0.1. The upper and bottom panels show the 
spectra produced in the jet and the jet + counter-jet, respectively. The 7-ray spectra above are shown for the minimum 
flux above 100 GeV, i.e. at the phase 0. counted from the location of the compact object behind the massive star (full 
curves), and for the maximum flux at the phase ~0.45 and 0.55 (dashed curves). The synchrotron spectrum produced 
by electrons in the jet are shown by the full curves below MeV energies. The spectrum from Cyg X-l reported by the 
COMPTEL telescope and the upper limit obtained based on the EGRET data (McConnell et al. 2000, 2002) are shown 
by thin full curves and marked by C and E, respectively. The sensitivities of the MAGIC and VERITAS telescopes 
are shown by double dot-dashed and dotted curves and marked by M and V, respectively (see e.g. Lorenz 2000). 



However, in the context of the recent discoveries of the 
TeV emission from massive binary systems (i.e. LS 5039 
and LS I +61 303), which have similar orbital parameters 
to Cyg X-l (i.e., a massive star and separation of the com- 
ponents), it is reasonable to predict the level and spectral 
features of the possible high energy 7-ray emission also 
from this object. Based on the IC e pair cascade model 
(Bednarek 2006b), considered here for Cyg X-l, we were 
able to predict the phase dependent TcV emission from 
LS I +61 303 (Albert ct al. 2006). In the case of Cyg X- 
1, we normalize the cascade 7-ray spectra expected from 
this object to the upper limit derived from the EGRET 
data. Note that, the MeV 7-ray emission from Cyg X- 
1 is ~ 1% of its total bolometric power (McConnell et 
al. 2002), and the EGRET upper limit in the GeV en- 
ergies is again on the level of 10% of the MeV emission. 
On the other hand, in the cascade model which we con- 
sider almost all energy from relativistic electrons in the jet 
is transfered to gamma-rays (and also secondary pairs), 
sec Bednarek (2006b). Since the jet power is expected to 
be comparable to the observed luminosity of the Cyg X- 
1, we conclude that the considered model is consistent if 



only a small part of jet energy is converted into relativistic 
electrons (i.e. ~ 0.1%). The applied normalization to the 
EGRET upper limit allows us to estimate the level of the 
TeV 7-ray flux and spectrum from Cyg X-l. In Figs. [6l 
we show the comparison of the results of calculations with 
the observational constraints at lower energies. The upper 
and the lower panels we show the cascade 7-ray spectra 
escaping toward the observer located at the inclination an- 
gle i = 30° for the case of their production in the jet and 
the jet+counter-jet (the sum of both), respectively. Since 
we have no detailed information about the spectral shape 
in the GeV energy range from the EGRET observations 
(but only the upper limit), electrons injected into the jet 
with different spectra and injection efficiencies along the 
jet have to be considered. At first we consider the model 
with very flat electron spectrum (differential spectral in- 
dex equal to -2) and the power transferred to electrons 
which changes with distance along the jet according to 
L e (z) oc z~ 2 (model A). In fact, such flat electron spec- 
trum seems to be not very realistic if we compare the 
spectral features of the source LS I +61 303 (the object 
of our previous modeling) and the spectrum of Cyg X- 
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1. The electron spectrum with the index 2 postulates the 
appearance of a completely new component in the high 
energy spectrum of Cyg X-l which does not lay on the 
extrapolation from the MeV energy range. The compar- 
ison of the calculated cascade spectra with the available 
data and the sensitivities of the present MAGIC telescope 
(M) and the future VERITAS telescope (V) is shown in 
Fig. [6^L,b. The dashed and full curves show the expected 
spectra when the compact object is in front of and behind 
the massive star HDE 226868. Clearly, for such model pa- 
rameters, Cyg X-l can be potentially detected at energies 
above ^100 GeV by the MAGIC telescope and above a 
few 100 GeV by the VERITAS telescope when the com- 
pact object is in front of the massive star. We also show 
the synchrotron spectrum which is produced by electrons 
accelerated in the jets assuming that the Doppler boosting 
is not important (bumps at lower energies marked by full 
curves). These spectra are at a much lower level than the 
observations of Cyg X-l in the MeV energy range and can- 
not be responsible for this emission. Note, however, that 
if significant Doppler boosting is discovered in Cyg X-l 
(the speed of the jet comparable to the velocity of light), 
the flux in the calculated synchrotron spectrum should 
go up with respect to the 7-ray spectra. This is due to 
the fact that the cascade 7-ray spectra are mainly formed 
outside the jet, i.e. inside the whole volume of the binary 
system where the plasma move with the wind velocity of 
the massive star. In such a case, the Doppler boosting 
can be only important for the synchrotron radiation but 
is clearly unimportant for cascade 7-rays. This MeV 7- 
ray emission can also originate in another scenario, e.g. 
in the hot corona of the accretion disk (e.g. Zdziarski & 
Gierlihski 2004, Chakrabarti & Mandal 2006) or in the 
inner part of the jet close to its base (a scenario not con- 
sidered in this paper) as a result of Comptonization of 
the disk radiation and possible IC cascading in the disk 
radiation, e.g. Aharonian et al. 1985, Coppi, Blandford & 
Rees 1993). 

The calculations with the more realistic model pa- 
rameters (see Table. 2) arc compared with the Cyg X-l 
data and sensitivities of the Chercnkov telescopes and are 
shown in Figs. [6] We show in Figs. Et the 7-ray spectra 
from the jet and in Figs. EH the jet + counter-jet, which 
are calculated for electrons with the differential spectral 
index oc E~ 2 - 6 , and the power in electrons dropping along 
the jet according to L c (z) cx z~ 2 (model B). For these 
parameters, there is no chance for detection of TeV 7- 
ray emission neither by the MAGIC telescope above ~100 
GeV nor by the Veritas telescope at ~1 TcV. Finally in 
Figs.[6j;f, the 7-ray spectra are shown for electrons injected 
with the steep spectrum (oc E~ 26 ) but with flat profile for 
the power transferred from the acceleration mechanism to 
electrons, L c (z) = const (model C). For this model, a 
large power is transferred to electrons at larger distances 
from the base of the jet where the acceleration of elec- 
trons is saturated at larger energies. Then, the cascade 
7-ray spectra at TeV energies are relatively flat, in spite 
of the strong 7-ray absorption effects. The 7-ray spectra 



Table 2. Parameetrs of the models considered in Fig. [SI 



Model ABC 

differential spectrum 

of electrons: N(E) oc E~ 2 g~ 2 ti g~ 2 fi 

injection number of electrons 

along the jet: L c (z) oc z~ 2 z~' 2 const. 



produced from the jet, after normalization to the EGRET 
upper limit in the GeV energy range, have chance to be 
detected by the VERITAS telescope at TeV energies but 
not by the MAGIC telescope above ~100 GeV (Figure 
6e). Note that the unboosted synchrotron spectra calcu- 
lated in model B and C are at the lower level than in the 
model A due to different normalization in the GeV energy 
range. 

6. Conclusions 

We constrain the parameters of the IC e ± pair cascade 
model for which we predict observable 7-ray fluxes from 
Cyg X-l at energies >100 GeV by the MAGIC and/or 
VERITAS telescopes. The main conclusions from the anal- 
ysis of this model for the 7-ray production in microquasars 
are reported in the following. 

The detection of the TeV 7-ray emission from Cyg X-l 
by Chercnkov telescopes can occur only under very spe- 
cial conditions for the acceleration of electrons in the jet as 
described in the model herein considered, and for an accel- 
eration efficiency parameter £, the magnetic field strength 
parameter 77, a flat spectrum of electrons, and a flat profile 
of the conversion of energy from the acceleration mecha- 
nism to electrons along the jet. These conditions are not 
very likely and do not fit nicely to 7-ray observations of 
Cyg X-l at the MeV- GeV energy range. 

Provided that such strict conditions for the accelera- 
tion of electrons are met (e.g. as in the considered model 
A), we have found that in the case of compact massive 
binary systems with circular orbits of the Cyg X-l type, 
the highest chances for detection of TeV 7-ray emission 
by the Chercnkov telescopes are for phases ~ 0.45 and 
~ 0.55, counted from the position of the compact object 
behind the massive star. The flux expected at the moment 
when the compact object is exactly in front of the mas- 
sive star is a factor of ^2 lower. The GeV and TeV 7-ray 
fluxes should be anticorrelated. The largest GeV flux is 
expected when the compact object is behind the massive 
star. The GeV 7-ray flux should vary by a factor of ~5 and 
the TeV flux by a factor of ~20 with the orbital period of 
Cyg X-1/HDE226868 system. 

We want to remind that the above conclusions con- 
cerning possible detection of the very high energy 7-ray 
emission from Cyg X-l are based on the assumption that 
the upper limit derived from the EGRET observations 
provides correct constraints on the multifrequency spec- 
trum of Cyg X-l, independent of the activity state of this 
binary system. In fact, this upper limit is consistent with 
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the extrapolation of the 7-ray spectrum reported by the 
COMPTEL telescope from the Cyg X-l region during the 
hard X-ray state of this source. Such flat MeV spectrum is 
not observed from Cyg X-l during its other activity states 
(e.g. McConnell et al. 2000,2002). 
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